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We evaluate heavy-quark (HQ) transport properties in a Quark-Gluon Plasma (QGP) within a 
Brueckner many-body scheme employing interaction potentials extracted from thermal lattice QCD. 
The in-medium T-matrices for elastic charm- and bottom-quark scattering off light quarks in the 
QGP are dominated by attractive meson and diquark channels which support resonance states up 
to temperatures of ~1.5 T c . The resulting drag coefficient increases with decreasing temperature, 
contrary to expectations based on perturbative QCD scattering. Employing relativistic Langevin 
simulations we compute HQ spectra and elliptic flow in y / ijvjv=200 GeV Au-Au collisions. A good 
agreement with electron decay data supports our nonperturbative computation of HQ diffusion, 
indicative for a strongly coupled QGP. 

PACS numbers: 25.75.-q 25.75.Dw 25.75.Nq 
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Experiments at the Relativistic Heavy-Ion Collider 
(RHIC) have shown that the matter produced in Au-Au 
collisions cannot be described by a weakly interacting gas 
of quarks and gluons, but rather consists of a strongly 
coupled Quark-Gluon Plasma (sQGP) with remarkably 
large opacity and low viscosity. The latter is required 
by hydrodynamic descriptions of the expanding fireball, 
implying rapid thermalization of the medium [J, Q . The 
understanding of these properties in terms of the under- 
lying interactions in the QGP, as governed by Quantum 
Chromodynamics (QCD), is a key theoretical objective. 
A valuable probe of the sQGP are heavy quarks (charm 
and bottom) which, due to their large mass, rriQ^>T c 
(T c ~180 MeV: critical temperature Q), are believed to 
be sensitive to the processes that establish and maintain 
thermalization of the medium, even at soft momentum 
scales. RHIC data for single-electron (e ) spectra as- 
sociated with semileptonic heavy-quark (HQ) decays in 
Au-Au collisions exhibit a surprisingly strong suppression 
and elliptic flow 0, [H, Q , indicating substantial collective 
behavior of charm quarks in the expanding fireball. Per- 
turbative QCD (pQCD) calculations, based on radiative 
energy loss, cannot explain these findings, even after in- 
clusion of elastic scattering 0, H[. Furthermore, it has 
been argued that the convergence of the perturbative se- 
ries for the HQ diffusion constant is rather poor Q , which 
calls for nonperturbative approaches. Effective models 
with strong HQ coupling in the QGP [H OH El El lead 
to significantly reduced thermal relaxation times com- 
pared to pQCD elastic scattering [Dm. resulting in better 
agreement [15|, EH with spectra [1, [H, [f| . 

In the present article, we perform a microscopic cal- 
culation of HQ diffusion in the QGP employing a non- 
perturbative T-matrix approach 17] with a driving ker- 
nel (potential) estimated from finite-temperature lattice 
QCD computations. We include a complete set of color 



channels for heavy-light quark interactions, as well as 
i=0,l partial waves together with HQ spin symmetry. 
This, in principle, provides an estimate of (elastic) trans- 
port coefficients without tunable parameters, albeit sig- 
nificant uncertainties remain in the definition of the po- 
tential. Within these uncertainties applications to HQ 
observables at RHIC support our approach. 

To evaluate in-medium properties of heavy quarks 
(Q=c,&) and heavy-light quark correlations, we employ 
a Brueckner-type many-body approach [ItJ ■ Our key as- 
sumptions are: (i) the main features of the elastic heavy- 
light quark interaction can be approximated by a static 
interaction potential, V(r) (to leading order in 1/mg, 
such an approach has been successfully applied for D 
meson spectra and decays in the vacuum [la . [l9j). and 
(ii) V(r) can be extracted from 1QCD simulations of the 
singlet-free energy F\ (r, T) [ 2g| f or a static QQ pair. As 
in previous works [l7l]2ll. |22|. |23||. we identify the poten- 
tial with the internal energy, U\, which is obtained by 
subtracting the entropy contribution from the free en- 
ergy, 



V 1 (r,T) = U 1 (r,T)-U 1 (oo,T), 



dFi 

Ui = Fi-T^. (1) 



A further subtraction is required to ensure the vanishing 
of the potential at large distance and thus the conver- 
gence of the T-matrix integral equation. In 1QCD sim- 
ulations the large distance limit of the internal energy, 
I7i (oo, T>T C ), is a decreasing function of the tempera- 
ture. It is tempting to associate this quantity with a 
selfenergy contribution to the HQ mass, m,Q(T)=mo + 
l7i(oo,T)/2 (mo: "bare" mass). However, around T c , 
C/i(oo, T) develops a rather pronounced maximum struc- 
ture rendering a mass interpretation problematic. Fur- 
thermore, little is known about the momentum depen- 
dence of this quantity. For simplicity, we assume constant 
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values for effective c- and 6-quark masses of m c =1.5 GeV 
and m;,=4.5 GeV (the difference to the current mass is 
mainly attributed to perturbative contributions). 

In addition to the color-singlet (meson) channel, we 
consider HQ interactions in the color-octet Qq, as well as 
in antitriplet and sextet Qq (diquark) channels. For the 
corresponding potentials we adopt Casimir scaling ac- 
cording to leading-order (LO) perturbation theory, Vs = 
— gVi, V3 = 5V1, Ve = — |Vi, which is also supporte d by 
1QCD calculations of the finite- T HQ free energy HHHf . 

The largest uncertainty in our calculations resides in 
the definition and extraction of the potential. While 
the identification with the internal energy (rather than 
the free energy) may be considered as an upper limit, 
the variations between different 1QCD calculations and 
pertinent parametrizations to numerically evaluate the 
entropy term in Eq. JT]), are appreciable. We have 
adopted 3 different potentials, based on parametrizations 
of Fi(r,T) in Refs. [H = [Wo], 0]=[SZ] and [IzjMMR] 
for quenched [H, 2-flavor [H and 3-flavor 1QCD [27j], re- 
spectively. The parametrizations [Wo] and [SZ] are simi- 
lar to a recent extraction [23| from 3-flavor 1QCD [28j]. 
The [MR] potential is deeper than the other two for 
T<1.6 T c (and consequently gives larger effects), but falls 
off faster above. The resulting transport coefficients vary 
by ^40%. More details will be elaborated in an extended 
paper 29]; in the following we restrict ourselves to the 
[Wo] potential. The lQCD-based potentials are imple- 
mented into a Brueckner many-body approach for heavy 
quarks, defined by a system of coupled Bethe-Salpeter 
(BS) and Dyson equations: 



T = K + J KGT , (2) 

E Q = E s + J TS\ S Q =S$ + S$H Q S Q , (3) 

with T: heavy-light quark T-matrix, K: interaction ker- 
nel, G: 2-particle propagator, S® q (S® ,q ): (free) single- 
particle propagators for heavy and light quarks, 
HQ selfenergy with contributions from thermal gluons 
(E ff ) and the T-matrix part from interactions with ther- 
mal light-anti-/quarks. Since we focus on a QGP at 
zero chemical potential (/Lt 9 =0), all quantities are quark- 
antiquark symmetric. To close the equations in the quark 
sector, one needs the corresponding system of equations 
for the light sector, which has been solved selfconsis- 
tently for E g and T qq in Ref. [l3]- Here, we augment 
these results by qq diquark interactions and implement 
the quark selfenergies in simplified form with constant 
real and imaginary parts as an effective quark mass, 
m q =0.25 GeV, and width, r g =0.2 GeV 17]]; variations 
in these parameters have little impact on the resulting 
quark selfenergies. The effects from heavy quarks in the 
heat bath can be safely neglected. 

To implement the (static) lQCD-based potential into 
our Brueckner approach we adopt the following approxi- 
mations. First, we use a standard reduction scheme (30j 
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sing, T=1.2T C 
sing, T=1.5 T c 
trip, T=1.2T C 
trip, T=1.5T C 
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FIG. 1: Imaginary part of the in-medium T-matrix for 
S— wave cq and cq scattering in color-singlet and antitriplet 
channels, respectively, at two different temperatures. The 
two-body threshold is at E t i 1Y =m c +m q — 1.75 GeV. 



to convert the 4-dimensional BS equation into a 3- 
D Lippmann-Schwinger (LS) equation (other reduc- 
tions [3l[ lead to very similar results |17|). thus neglect- 
ing virtual particle-antiparticle loops but keeping rela- 
tivistic kinematics of the quark propagators. This, in 
turn, enables to identify the reduced kernel K with the 
potentials V 1 3 e 8 constructed above, representing a lad- 
der approximation to the T-matrix. As in Refs. [17l. [2]|. 
we correct the static potentials for a relativistic Breit 
(current-current) interaction. Azimuthal symmetry and 
a partial-wave expansion lead to a 1-D LS equation in 
each color (a) and angular-momentum (I) channel (E: 
center-of-mass energy of the Qq system) , 



T a , l (E;q',q) = V a ,i(q',q) + - J dk k 2 V a ,i(q',k) 
xG Qq (E; k)T a , t (E; k, q)[l - - . 



/(to) denote Fermi-Dirac distributions and u>\ — (mf + 
fc 2 ) 1 / 2 quasiparticle dispersion laws. We include both 5*- 
(1=0) and P-wave {1=1) channels. The 2-particle propa- 
gator in the Thompson scheme [3(| reads 



G(£7;fc) = (l/4)/[£;-( W «+iE«) 



(u£+i£?)]. (5) 



Results for the in-medium S-wave T-matrix are illus- 
trated in Fig. [T] for c-quark scattering. The attractive 
color-singlet and -antitriplet channels are the dominant 
contributions, supporting resonance structures up to 
temperatures of ~1.7T C and ~1.4T C , respectively. Both 
the repulsive color channels, as well as P- waves, lead to 
much smaller T-matrices. However, due to larger de- 
generacies their contribution to the HQ selfenergies and 
transport coefficients is not negligible. 

Next we compute the HQ selfenergies and transport 
coefficients generated by the nonperturbative heavy- light 
T-matrices. The HQ selfenergy represented by the sec- 
ond term on the right-hand-side of Eq. ([3]) is calculated 
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FIG. 2: Real (upper lines) and imaginary (lower lines) parts 
of the on-shell c-quark selfenergy versus 3-momentum at tem- 
peratures T=l.l, 1.4 and 1.8 T c . 
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FIG. 3: Friction coefficients of charm quarks vs. 3-momentum 
for different temperatures in a QGP, based on our nonpertur- 
bative T-matrix approach (upper curves at p=0) and on LO 
pQCD scattering with a s =0.4 (lower curves). 



within the imaginary-time formalism as 
dsida 



£?C*»;p) 



G 



dV 



(2nf 



(-T) 



T a (z v + 2v;p, p')£)q(^',p') 



(6) 



(z u =\tt(2v+1)T: fermionic Matsubara frequencies). As 
implicit in our T-matrix (potential) we assume spin 
and light-flavor degeneracy of the heavy-light interaction 
(in line with the free P/-meson spectrum [32|), yield- 
ing d S i=A{l2)N f for S'(T)-waves (iV/=2.5 to account 
for the smaller strange-quark density). The resulting c- 
quark selfenergies (Fig. [2|) translate into large in-medium 
widths of around T c =-2 ImE c ~200 MeV (consistent with 
our input parameters). The dominant meson and diquark 
contributions are about equal, while the P-wave amounts 
to ~40%. The nonperturbative real parts are small. 
We evaluate HQ diffusion in the QGP within a Fokker- 



Planck equation for the distribution function, /q [lij 

d 2 



dpidp 



■(By fa) , 



(7) 



with momentum and temperature dependent friction (7) 
and diffusion (P^) coefficients. They are calculated as 
in Ref. [lfj using the appropriate relation between the 
invariant amplitude Ai and our in-medium T-matrix, 



64tt 



I 2 1 2 \2 / 2 1 2\1 

(s -m + m Q ) {s-m Q + m ) 



xN f Vda(|T a , /=0 (s)| 2 + 3|T a , i=1 (s) cos(# cm )| 2 ) . 



(8) 



The non-perturbative thermal relaxation rates reach up 
to 7~l/(7fm/c) at low momenta close to T c , a factor 
of ~4 larger than elastic pQCD scattering (but compa- 
rable to the resonance model of Ref. [lfj), cf. Fig. [3] 
In contrast to other calculations available thus far, the 
thermalization rate decreases with temperature, due to 
the dissolving resonances induced by the screening in the 
lQCD-based potentials. The increase in quark density 
is overcompensated by the loss of interaction strength. 
This has important consequences for HQ observables at 
RHIC, as discussed below. The cos 2 (# cm ) factor in \M\ 2 
reduces the P-wave contribution to 7 to ^20% of the 
S'-wave. Combining T-matrix and pQCD contributions, 
the spatial HQ-diffusion constant at p—0 amounts to 
DHQ=T/(m c 7)~5/27rT at T=200 MeV, a factor of -4 
smaller than in pQCD, thus corroborating the notion of 
a strongly coupled QGP at temperatures up to 1.5-2 T c . 

The above calculated transport coefficients (from the 
T-matrix plus LO pQCD scattering off gluons) are im- 
plemented into Langevin simulations of c and b quarks in 
Au-Au collisions at RHIC using an isentropically expand- 
ing QGP fireball [l6(. The latter has been constructed 
to resemble hydrodynamic models, with a bulk elliptic 
flow of ~5.5% and initial temperature of T o =340 MeV 
in semicentral collisions (when using LO pQCD interac- 
tions, the resulting HQ spectra are in good agreement 
with Langevin simulations in an explicit hydrodynamic 
background [TlT|). To compare to experimental spec- 
tra, c and b quarks are hadronized at the end of the QGP- 
hadron gas mixed phase within the coalescence model of 
Ref. [331 ] (supplemented by 5-function fragmentation and 
D- and B- meson decays). Note that the resonance cor- 
relations in the T-matrix naturally merge into a quark- 
coalescence description toward T c [34j. Our calculations 
for the e^- nuclear modification factor, Paa (defined as 
the ratio of the spectrum in Au-Au collisions to the one 
in p-p scaled by the number of binary N-N collisions), 
and elliptic flow coefficient, V2, show fair agreement with 
recent RHIC data [U @ , cf . Fig. 0J Inspection of the time 
evolution of the c-quark distribution reveals that the sup- 
pression in the pr spectra (i.e., Raa) is mostly built up 
in the early stages of the QGP, while most of the vi is 
developed at temperatures close to T c . This feature is 
amplified by the temperature dependence of the trans- 
port coefficients in the T-matrix approach, and seems to 
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FIG. 4: Our results for the nuclear modification factor (up- 
per panel) and elliptic flow (lower panel) of single electrons 
with/without (solid/dashed lines) quark coalescence in Au-Au 
collisions compared to RHIC data 5, 6j. The estimated the- 
oretical uncertainty due to different extractions of the 1QCD 
based heavy-light quark potentials is up to ~30% [291 ]. 



be favored by the data (coalescence with light quarks 
further contributes to the increase in both t>2 and Raa)- 



In summary, we have calculated HQ selfenergies and 
transport coefficients within a T-matrix approach for 
heavy-light quark interactions in the QGP using two- 
body potentials estimated from lattice QCD. HQ scatter- 
ing turns out to be dominated by "prehadronic" mesonic 
and diquark channels which increase in strength when ap- 
proaching T c . These correlations substantially accelerate 
thermal relaxation times compared to pQCD and provide 
for a natural onset of the hadronization process (in such 
a scenario, nonperturbative HQ interactions with gluons 
are less relevant). When implemented into Langevin sim- 
ulations at RHIC, reasonable agreement with the sup- 
pression and elliptic flow of spectra from HQ decays 
emerges. This is rather remarkable in view of the largely 
parameter-free calculation of the transport coefficients. 
Future work should aim at scrutinizing the uncertainties 
inherent in the potential approach at finite temperature 
and in the extraction of the potential from lattice QCD. 
Further insights could be obtained from direct 1QCD 
computations of heavy-light quark correlation functions 
in the QGP. In addition, elastic HQ interactions, which 
parametrically dominate at low pt, should be supple- 
mented by radiative energy loss [131 ] which takes over at 
high pt- Our present analysis suggests that a small HQ 
diffusion coefficient arises from a nonperturbative inter- 
action strength in a strongly coupled QGP. 
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